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INTRODUCTION 


The  primary  factors  affecting  infrared  electro-optical  surveillance,  guidance  and 
weapon  systems  in  the  marine  environment  are  atmospheric  water  vapor  and  aerosols 
which  absorb  and  scatter  the  radiation.  In  the  absence  of  real-time  measurements,  we 
must  presently  rely  on  the  LOWTRAN  6  atmospheric  propagation  code1  to  predict 
infrared  transmission  losses  and  sky  backgrounds,  using  as  inputs  measured 
meteorological  parameters.  The  effects  of  water  vapor  absorptions  are  adequately 
handled  by  LOWTRAN  6.  and  selectable  size  distribution  models  are  available  for 
calculating  the  aerosols'  absorption  and  scattering  properties.  One  of  these  aerosol 
models  (Navy  Maritime  Model),  which  is  applicable  to  ocean  atmospheres,  was 
developed  by  Gathman2  at  the  Naval  Research  Laboratory,  utilizing  a  large  data  set  of 


size  distributions  and  meteorological  parameters  measured  near  the  surface  in  a  variety 
of  marine  environments.  The  particle  size  distribution  model  (at  radius  r)  is  the  sum 
of  three  log-normal  distributions  given  by 

n(r)  =2  Ai  exP  ['(*?/)  ]'  cm  /*"»  ' 

i=l  ' 

(1) 

where 

Aj  =  2000  (AM)2 

(2a) 

A2  =  5.866  (v-2.2) 

(2b) 

A3  =  0.01527  (vc-2.2)  . 

(2c) 

Component  one  represents  the  contribution  by  continental  aerosols.  (AM)  is  an  air 
mass  parameter  which  varies  between  integer  values  of  1.0  for  open  ocean  to  10  for 
coastal  areas  given  by 


(AM)  =  Rn/4  ♦  1,  (3) 

where  Rn  is  the  measurement  of  atmospheric  radon  content  expressed  in  picocuries 
per  cubic  meter  (pCi/m3).  Components  two  and  three  represent  equilibrium  sea  spray 
particles  generated  by  the  24-hour  averaged  (v)  and  current  (vc)  surface  wind  speeds 
in  meters  per  second.  In  Eq.  (2b).  if  5.866(v-2.2)  <  0.5.  then  A2  =  0.5:  and  in  Eq. 
(2c).  if  0.01527(vc-2.2)  <  1.5  x  10‘5.  then  A3  =  1.5  x  10  5. 

In  Eq.  (1).  n  .  the  modal  radius  for  °ach  comnonent.  is  allowed  to  grow  with 
relative  humidity  (RH)  according  to 

'  ■  M-*y]1/3  •  «> 


The  contribution  to  the  total  extinction  or  absorption  by  each  aerosol  component  can 
be  written  as 

/».,.( M-,  =  Cj  I  .xp  .  (5) 

where  Ci  =  The  factor  f-1  in  the  expression  for  C)  insures  a  constant 

total  number  of  particles  as  the  relative  humidity  increases.  Q*,a(X.r.m)  is  the  cross 
section  for  either  extinction  or  absorption  normalized  to  the  spherical-particle 
geometrical  cross  section,  and  m  is  the  complex  refractive  index,  which  is  allowed  to 
change  from  that  of  dry  sea  salt  as  the  particle  deliquesces  with  increasing  humidity. 
LOWTRAN  6  provides  precalculated  values  in  tabular  form  of  the  parameter  /J*,*(X)j/C| 
at  discrete  wavelengths  and  four  relative  humidities  (50,  85.  90  and  99%).  from  which 
the  average  extinction  or  absorption  coefficient  for  a  specific  wavelength  band  and 
relative  humidity  can  be  readily  determined  by  interpolation  between  the  stored  values. 
When  an  observed  surface  visual  range  (visibility)  is  available  as  an  input  to  the 
model,  the  amplitudes  of  the  three  components  will  be  adjusted  so  that  the  calculated 
visual  range  at  a  wavelength  of  0.55  Jim  is  the  same  as  the  observed  value. 

The  accuracy  to  which  this  model  can  predict  infrared  extinction  coefficients  has 
been  tested  only  against  a  limited  set  of  surface  transmissometer  and  meteorological 
measurements  at  San  Nicolas  Island3.  Good  correlations  between  calculated  and 
measured  extinctions  for  wavelengths  of  1.06  /im  and  3.6  /<m  were  obtained.  At  10.5 
ft m  the  agreement  was  less,  with  the  calculated  extinctions  being  20  to  40%  greater 
than  those  measured  by  the  transmissometer.  These  correlations,  however,  were 
sensitive  to  the  selection  of  the  air  mass  factor  and  whether  or  not  the  visibility  was 
used  as  an  input.  An  alternative  approach  is  to  test  the  model's  utility  to  predict 
the  infrared  radiance  of  the  sky.  It  is  well  known4  that  the  absorption  (and 
emissivity)  of  the  atmosphere  for  the  8-  to  12  pm  wavelength  band  depends  on  the 
optical  path  length  such  that  the  effective  blackbody  temperature  of  the  sky  will 
increase  with  the  zenith  angle.  Near  the  horizon,  the  sky  temperature  will  equal  the 
ambient  air  temperature  unless  aerosols,  which  scatter  the  radition,  are  present.  The 
effects  of  aerosols  for  this  wavelength  band,  however,  are  noticeable  only  at  zenith 
angles  greater  than  about  85  deg  (in  cloud-free  skies).  In  this  paper,  we  examine  the 
utility  of  the  aerosol  model  (with  the  LOWTRAN  6  radiance  algorithm)  to  predict 
infrared  (8  -  12  /<m)  sky  radiances  which  were  measured  close  to  the  horizon 
simultaneously  with  radiosonde  measurements  of  meteorological  parameters.  For  this  to 
be  a  valid  approach,  we  must  rely  upon  the  accuracies  of  the  measurements  and  the 
LOWTRAN  6  radiance  algorithm. 

MEASUREMENTS 

The  infrared  (8  -  12  /im)  sky  radiances  for  these  investigations  were  obtained  on 
16  April  1986.  The  measurements  were  made  with  a  calibrated  thermal  imaging 
system  (AGA  THERMOVISION.  Model  780)  using  a  2.9?  field-of-view  (FOV)  Tens 
with  an  instantaneous  field-of-view  ( I F 0 V)  of  0.9  mrad.  The  response  of  each 
wavelength  band  is  determined  by  placing  a  blackbody  of  known  temperature  (  ±  0.1°  C 
for  temperatures  <5C f  C)  close  to  the  aperture  of  the  lens.  The  digitized  video  signal 
transfer  function  of  the  system  then  allows  the  blackbody  temperature  to  be 
reproduced  to  within  0.2°  C.  The  video  output  of  the  scanner  is  digitized  and 
processed  on  a  microcomputer  to  allow  the  temperature  of  selected  pixels  of  the  scene 


to  be  displayed.  For  these  measurements  the  scanner  was  directed  due  west  over  the 
ocean  from  an  altitude  of  33  m  such  that  approximately  T  of  the  FOV  was  above 
the  horizon.  During  the  recording  period,  radiosondes  were  launched  from  a  ship  (USS 
Point  Loma  (AGDS-2))  5  km  off  the  coast  of  Pt.  Loma.  San  Diego.  CA.  The 
radiosonde  system  employed  was  the  VAISALA  model  RS80.  The  measured 
temperature  and  relative  humidity  variations  with  altitude  for  three  periods  0845.  1245 
and  1645  PST  16  April)  are  graphically  shown  in  Fig.  1  and  tabulated  with  the 
pressure  variations  in  Table  1.  During  the  first  launch,  broken  stratus  clouds  were 
present  near  an  elevation  of  900  m.  During  the  subsequent  launches,  the  clouds 
persisted,  but  the  coverage  was  scattered.  Surface  wind  speeds  and  directions  were 
recorded  continuously  on  shore  at  the  sensor  site  and  periodically  aboard  the  ship. 
Northwesterly  winds  (31CP  ±  1CP  )  had  persisted  for  24  hours  prior  to  and  during  the 
measurements  with  varying  speeds  as  shown  in  Fig.  2.  At  0800  PST  on  16  April 
the  wind  speed  had  increased  from  approximately  3  m/s  to  values  between  9  m/s 
and  12  m/s.  The  24-hour  average  and  current  wind  speeds  coinciding  with  the  times 
of  the  radiosonde  launches  are  tabulated  in  the  figure.  Measurements  of  atmospheric 
radon  were  also  made  aboard  the  USS  Point  Loma  to  aid  in  determining  the  air  mass 
characteristics.  The  radon  counts  measured  as  a  function  of  time  are  shown  in  Fig. 
3  and  indicate  the  air  mass  was  primarily  of  maritime  origin  (<4  pCi/m!)  throughout 
the  measurement  period.  The  increased  radon  counts  near  0400  PST  on  15  April 
coincide  with  the  in-port  time  of  the  ship. 


COMPARISON  OF  MEASUREMENTS  AND  CALCULATIONS 


The  LOWTRAN  6  radiance  algorithm  assumes  the  atmosphere  to  be  composed  of 
a  number  (nmax  =  33)  of  isothermal  layers  characterized  by  a  temperature  Tj  and 
spectral  transmittance  r(X.  i.  8)  along  the  optical  path  traversing  the  ith  layer  at  angle 
8  From  Kirchoff's  law.  the  spectral  radiance  of  the  i'h  layer  is 


N(A,i,B)  =  1-r  (A,a,8) 


B(X,T.) 


(6) 


where  ra(X.i.8)  is  the  absorption  transmittance  and  B(X.Tj)  is  Planck  s  blackbody 
radiation  formula  Then  the  spectral  radiance  reaching  the  ground  through  the 
intervening  atmosphere 


i-1 


N(X,i,B)  n  T(X, i ,8)  =  [l  -  ra{A,i,B) 
j=l 


1  r'  -1 

n 

j=i 


T  B( X ,  T  - ) 


(7) 


Summing  up  the  contribution  from  all  layers,  the  spectral  radiance  reaching  the  ground 
is 
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In  t  h  t  s  equation  scattering  is  considered  only  as  a  loss  mechanism  through  the 
extinction  transmittanr  e  term  1 1  'T ( X  i.H)  and  is  not  included  as  a  source  ot  radiation 
It  lias  In  cn  proposed  by  Ben  Shalom  i  t  al  1  that  the  LOWTRAN  algorithm  was 
debt  lent  ie  multiple  scattering  effects  over  the  long  propagation  paths  affecting  the 
sky  radium  e  were  not  properly  addressed  They  proposed  a  modification  to  LOWTRAN 
to  include  scattering  as  a  source  of  radiation  by  replacing  the  absorption  transmittance 
m  Eq  (8)  hv  the  extinction  transmittance  However,  utilizing  data  similar  to  that 
lu  icin  Hughes  it  .i  I  *  have  shown  that  the  proposed  conservative  scattering 
modifications  to  LOWTRAN  grossly  overestimate  the  horizon  sky  radiances  when 
aerosols  .ire  present  and  tfiat  multiple  scattering  effects  are  negligible,  at  least  for  the 
wavelength  band  and  atmospheric  conditions  considered 

In  fig  4  the  sky  radiances  measured  at  the  optical  horizon  (zenith  angle.  8  = 
90  17°  I  ire  compared  to  those  calculated  using  the  unmodified  LOWTRAN  6  code 
with  the  measured  meteorological  data  We  have  chosen  to  address  only  the  sky 
radiance  at  t  ho  optical  horizon  because  of  the  possible  contamination  of  the 
measurements  by  the  scattered  stratus  clouds  (It  can  be  shown'  that  8  —  12  /cm 
radiances  at  the  optical  horizon  are  insensitive  to  cloud  emissions  because  of  the  low 
atmospheric  transmittanc es  over  the  contributing  optical  path  lengths.)  The  clear-air 
radiance  calculations  were  made  using  plus  and  minus  uncertainties  (0  5°  C  in 
temperature  and  5%  in  relative  humidity)  as  shown  In  each  case,  the  clear  air 
calculations  are  greater  than  the  mea  .urements.  indicating  a  small  presence  of 
aerosols  (These  radiance  differences  correspond  to  equivalent  blackbody  temperature 
differences  of  2  to  3°  C  )  The  calculations  made  with  the  Navy  Maritime  Aerosol 
Model  were  for  an  air  mass  factor  of  unity  for  maritime  air  (as  indicated  by  the 
radon  measurements)  and  the  24  hour  average  and  current  surface  wind  speeds  as 
listed  in  Fig  2  In  the  first  case  there  is  good  agreement  between  the  measured  and 
calculated  radiances  By  adjusting  the  surface  visibility  input  to  130  km  (as  compared 
to  the  default  value  of  96  8  km)  the  calculated  radiance  can  be  made  to  coincide  with 
the  measured  value  For  the  second  and  third  time  periods,  the  calculations  differ 
greatly  from  tin  measurements  by  equivalent  blackbody  temperatures  of  approximately 
1 9°  C  and  20°  C.  respectively  The  calculations  can  be  made  to  agree  with  the 

measurements  by  adjusting  the  default  visibilities  of  31  km  and  35  km  to  values  of 

180  km  and  210  km  respectively  These  visibilities  are  excessive,  based  on  visual 

observations  ol  coastal  islands  at  the  time  of  the  measurements  Los  Coronados 
Islands  (-30  km  distant)  were  clearly  seen  However.  San  Clemente  Island  with  a  peak 
•  legation  of  ~600  m  was  not  visible  from  the  upper  decks  of  the  USS  Point  Loma  at 

a  distance  of  75  km  In  the  8-12  ^m  band,  the  horizon  radiance  is  affected  mainly 

by  the  aerosols  with  radii  greater  than  1  /itn.  which  scatter  the  radiation.  The 
comparison  discrepancies  in  Fig  4  most  likely  stem  from  the  current  wind  speed 
component  which  determines  the  number  of  particles  greater  than  1  /jm.  This 
wavelength  hand  is  less  likely  to  be  influenced  by  the  24  hour  average  wind-speed 

component  which  mainly  generates  particles  in  the  0  1  to  1  -/im  radius  interval.  In 

Fig  5  the  relative  sensitivity  of  the  radiance  calculations  to  the  wind  speed  factors  is 

demonstrated  by  means  of  the  1645  PST  data  set  The  radiance  calculations  are 

insensitive  tc  ?4  hour  wind  speeds  varying  between  2  2  m/s  and  10  m/s,  but  are 
extremely  serritive  to  the  current  wind  speed  If  the  multiplying  constant  in  Eq.  (2c). 
(k  0  01527  s'rn)  and  a  24  hour  average  wind  speed  of  6.0  m/s  are  maintained, 
the  current  wind  speed  must  be*  reduced  to  2  6  m/s  to  obtain  agreement  between 
tabulated  and  measured  radiant  is  If  the  measured  value  of  current  wind  speed  is  to 
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Figure  4  Comparison  of  infrared  sky  radiances  measured  at  the  optical  horizon  with  those  calculated 
using  LOWTRAN  6. 


Figure  5  Sensitivity  of  the  optical  horizon  sky  radiances  calculated  with  LOWTRAN  6  using 
the  Navy  Maritime  Aerosol  Model  with  differing  surface  wind-speed  parameters. 
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be  maintained  in  the  LOWTRAN  calculations,  the  multiplying  constant,  k.  must  be 
reduced  by  a  factor  near  20  to  0.0007  s/m.  The  size  distributions  associated  with  the 
different  values  of  k  are  shown  in  Fig.  6.  For  particle  radii  less  than  1  pm.  the  size 
distributions  are  unaffected  by  changes  in  k.  At  larger  radii,  the  size  distribution  for 
k  =  0.0007  s/m  are  smaller  by  an  order  of  magnitude.  This  reduction  in  the  number 
of  larger  particles  accounts  for  the  increase  in  the  calculated  radiance,  i.e..  the 
scattering  losses  are  reduced. 

Recently,  published  results  by  Dee  Leeuw8  of  large  particle  size  distributions  (r  > 
5  /*m)  in  the  North  Atlantic  provide  another  method  of  evaluating  the  current  wind- 
speed  component  in  the  present  model.  In  that  work,  size  distributions  were 
measured  with  an  impactor  at  different  heights  (0.2  to  11  m)  above  the  sea  surface. 
Measured  size  distributions  (normalized  to  a  relative  humidity  of  80%  according  to  the 
formulas  of  Fitzgerald9)  were  graphically  presented  for  an  altitude  of  11  m  as  a 
function  of  surface  wind  speed.  Here,  size  distributions  were  calculated  with  the 
model,  using  the  measured  surface  meteorological  parameters  and  the  k-factor  in  the 
original  model  was  adjusted  to  obtain  agreement  with  those  presented  by  Dee  Leeuw 
(after  adjustment  to  the  average  measured  relative  humidity  in  the  first  two 
radiosonde  levels).  An  example  of  the  comparisons  is  shown  in  Fig.  7  for  the  1645 
PST  set  of  data.  Excellent  agreement  between  the  adjusted  and  measured  size 
distribution  is  obtained  for  a  k-  factor  of  0.00109  s/m.  This  value  is  approximately 
36%  higher  than  that  determined  from  the  radiance  measurements.  The  difference  may 
reflect  the  assumption  in  the  present  model  that  the  number  of  surface-generated 
particles  remains  constant  up  to  an  altitude  of  2  km.  where  the  LOWTRAN  6 
calculations  default  to  the  Tropospheric  Aerosol  Model.  However,  it  can  be  shown 
that  the  radiance  calculations  at  the  optical  horizon  are  affected  less  than  2%  by 
including  only  the  lowest  two  levels  of  the  radiosonde  profile.  The  k-factors 
determined  by  both  techniques  are  shown  in  Fig.  8  at  the  measured  wind  speeds. 
Within  the  measurement  accuracies  of  both  techniques  and  those  to  which  the  size 
distributions  could  be  scaled  from  the  graph  in  Dee  Leeuw’s  paper,  the  k-factors  can 
be  considered  to  be  in  reasonable  agreement. 


DISCUSSION 

It  is  interesting  to  notice  in  Fig.  6  the  suggestion  of  a  linear  dependency  of  the 
factor  k  on  the  current  wind  speed.  However,  because  of  the  small  data  sample,  no 
quantitative  conclusions  can  be  made  in  this  regard. 

A  joint  effort  by  the  Naval  Research  Laboratory,  the  Naval  Postgraduate  School, 
and  the  Naval  Ocean  Systems  Center  is  presently  underway  to  develop  a  Navy  Ocean 
Vertical  Aerosol  Model  (NOVAM)  for  inclusion  into  a  future  version  of  LOWTRAN. 
Using  the  current  LOWTRAN  6  Navy  Maritime  Aerosol  Model  as  the  surface  kernel, 
this  new  model  is  intended  to  greatly  reduce  the  third  component's  variation  with 
altitude.  The  results  of  this  study,  however,  have  demonstrated  that  for  moderate 
wind-speed  conditions,  the  current  wind-speed  component  in  the  kernel  model  may  be 
factors  near  20  too  large.  Therefore,  a  careful  re-examination  should  be  given  to  the 
constants  of  the  present  model  before  inclusion  to  NOVAM. 
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Figure  6.  Examples  of  aerosol  size  distributions  calculated  with  different  values 
of  the  multiplying  factor  k. 


CURRENT  WIND  SPEED  (m/s) 


Figure  8.  Values  of  the  current  wind-speed  multiplying  factor  k  required 
to  obtain  agreement  between  measured  and  calculated  radiances  and 
particle  size  distributions 
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